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• ^It  appear*  that  nonlinear  effect*  are  aore  severe  at  a 
speed- to-length  ratio  of  6 than  of,  sap,  4 or  las*,  principally 
because  of  the  reduction  of  the  doping  ratio  of  the  boat  with 
increasing  spaed,  and  the  consequent  increase  in  notions  in  the 
vicinity  of  the  resonant  encounter  frequency.  However,  it  IS 
concluded  that  the  linear  theory  in  its  present  fom  can  serve 
as  a useful  design  tool,  especially  since  it  has  already  been 
shown  that  the  effect  of  hull  paraneters  on  performance  as 
obtained  from  data  in  the  linear  range  is  valid  for  operation 
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ABSTRACT 

A theoretical  method  is  derived  for  predicting  the  linearized 
response  characteristics  of  constant  deadrise,  high-speed  planing 
boats  in  head  and  following  waves.  Comparisons  of  the  theoretical 
predictions  of  the  pitch  and  heave  response-amplitude  operators 
and  phase  angles  with  existing  experimental  data  show  reasonably 
good  agreement  for  a wide  variety  of  conditions  of  interest. 

It  appears  that  nonlinear  effects  are  more  severe  at  a speed- 
to-length  ratio  of  6 than  of,  say,  4 or  less,  principally  because 
of  the  reduction  of  the  damping  ratio  of  the  boat  with  increasing 
speed,  and  the  consequent  increase  in  motions  in  the  vicinity  of 
the  resonant  encounter  frequency.  However,  it  is  concluded  that 
the  linear  theory  in  its  present  form  can  serve  as  a useful 
design  tool,  especially  since  it  has  already  been  shown  that  the 
effect  of  hull  parameters  on  performance  as  obtained  from  data  in 
the  linear  range  is  valid  for  operation  in  realistic  seas. 
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This  investigation  was  authorized  and  funded  by  the  Naval  Sea  Systems 
Command  (SEA  035)  under  the  General  Hydrodynamics  Research  Program, 
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INTRODUCTION 

The  history  of  the  development  of  planing  hulls  has  depended  almost 
exclusively  in  the  past  on  acquisition  and  analysis  of  full-scale  and 
model  experimental  data.  This  has  been  especially  true  in  the  area  of  sea- 
keeping where  until  only  a few  years  ago,  practically  no  experimental  data 
were  available.  Following  a study  by  Savitsky,^  in  which  he  noted  that  the 
lack  of  such  data  led  to  continuous  controversy  among  boat  designers  as  to 

the  makings  of  a good  rou^h  water  boat,  extensive  experiments  were  carried 

2 

out  at  the  Davidson  Laboratory.  In  1969  Fridsma  carried  out  experiments 
on  a series  of  constant-deadrise  models  in  smooth  water  and  regular  head 
waves  to  define  the  effects  of  deadrise,  trim,  loading,  speed,  and  length- 
to-beam  ratio  as  well  as  wave  proportions  on  added  resistance,  heave  and 
pitch  motions,  and  impact  accelerations.  Despite  the  fact  that  he 


^Savitsky,  D. , "On  the  Seakeeping  of  Planing  Hulls,"  Marine  Technology 
(Apr  1968).  A complete  listing  of  references  is  given  on  page  82. 

2 

Fridsma,  G. , "A  Systematic  Study  of  the  Rough-Water  Performance  of 
Planing  Boats,"  Davidson  Laboratory,  Stevens  Institute  of  Technology 
Report  R-1275  (Nov  1969). 
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observed  a significant  nonlinear  dependence  of  the  motions  on  wave  height, 

he  concluded  that  his  findings,  which  were  predominantly  obtained  from 

measurements  in  the  linear  range,  were  valid  and  extremely  important. 

This  was  further  confirmed  by  later  experiments  made  with  more  realistic 

3 

irregular  head  waves,  when  he  observed  that  the  results  generally 

correlated  well  with  his  regular  wave  findings. 

The  relative  success  of  a recently  developed  theoretical  method  for 

4 

predicting  porpoising  of  planing  hulls  gave  encouragement  to  the  belief 
that  a reasonably  good  dynamic  model  for  the  planing  hull  was  in  hand  for 
the  first  time.  This  consisted  of  a set  of  three  linear  equations  in  the 
surge,  heave,  and  pitch  motions  of  the  boat  in  calm  water.  The  solution 
of  these  homogeneous  equations  for  the  stability  roots  led  to  predictions 
of  porpoising  trim  angles  which  were  on  the  whole  in  good  agreement  with 
experimental  data.  Expressions  for  the  linearized  wave  forces  and 
moments  due  to  regular  waves  were  consequently  derived  in  the  present  report 
and  were  used  with  the  previously  described  equations  for  predicting 
motions  in  waves. 

Although  the  theory  is  linear  and  has  several  approximations,  it  is 

felt  that  its  basic  structure  is  correct  in  view  of  the  good  overall 

2 

agreement  with  the  extensive  model  data  obtained  by  Fridsma  in  regular 
waves.  Furthermore,  as  noted  previously,  since  results  obtained  in  the 
linear  range  are  valid,  concerning  the  effect  of  hull  parameters  on 
performance,  it  is  felt  that  the  theory  will  serve  as  a useful  design  tool 
in  its  present  form. 
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EQUATIONS  OF  MOTION 

The  linearized  equations  for  the  longitudinal  motions  of  surge,  heave 
and  pitch  in  regular  head  and  following  waves  are  derived  in  Appendix  A. 
Although  the  derivation  is  for  all  three  degrees  of  freedom,  the  numerical 


3 

Fridsma,  G.,  A Systematic  Study  of  the  Rough-Water  Performance  of 
Planing  Boats,  Irregular  Waves — Part  2,"  Davidson  Laboratory,  Stevens 
Institute  of  Technology  Report  R-1495  (Mar  1971). 

4 

Martin,  M. , "Theoretical  Prediction  of  Porpoising  Instability  of  High- 
Speed  Planing  Boats,"  DTNSRDC  Report  76-0068  (Apr  1976). 
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calculations  in  this  report  have  been  obtained  from  the  solutions  to  the 

coupled  pitch  and  heave  equations  only.  This  was  felt  to  be  justified  on 

2 3 

the  basis  of  the  experimental  findings  of  Fridsma.  * He  found  from  visual 
observations  and  on  examination  of  the  time  history  record  of  model  boat 
motions  in  regular  head  waves  that  little  surging  motion  took  place, 
particularly  at  speed-to-length  ratios  of  4 and  6.  To  verify  this  observa- 
tion at  a speed-to-length  ratio  of  2,  comparative  tests  were  made  for  both 
"constant  speed"  and  "constant  thrust"  conditions.  He  found  little 
difference  in  the  pitch  and  heave  motions.  This  result  is  not  altogether 
surprising  in  view  of  the  fact  that,  according  to  Equation  (121)  of 
Appendix  A,  the  theoretical  estimates  of  the  surge  component  of  the  wave 
excitation  force  are  essentially  given  by  the  product  of  the  heave 
component  with  the  tangent  of  the  steady  state,  running-trim  angle.  Since 
the  trim  angles  of  his  experiments  were  between  4 and  6 degrees,  the  wave 
excitation  force  in  surge  would  be  expected  to  be  much  less  than  in  heave. 
Accordingly,  the  surge  degree  of  freedom  is  omitted  from  the  equations  of 
motion  given  as  follows.  These  equations,  which  are  in  nondimens ional 
form,  are  for  a planing  boat  moving  with  constant  forward  speed  U and  are 
written  with  respect  to  its  body  coordinate  system  at  its  steady  calm 
water  running  trim;  i.e.,  no  motion  except  steady  forward  speed.  The 
origin  of  the  coordinate  system  is  located  at  the  boat  center  of  gravity. 
The  axes  Ox,  Oy,  and  Os  are,  respectively,  forward,  starboard,  and 
vertically  down.  From  Equations  (127)  and  (128)  of  Appendix  A,  we  have 


(m'-Z..')  z’  - Z.'z’  - Z z'  - Zv'Q'  - Za'0'  - Z ' 0 = (Z)  ’ e C 

Z Z Z 0 0 0 w 

-iw  't' 

- M..'z'  - M. ’z'  - M 'z'  + (I  '-Ms')  0'  - Ma'0'  - Moe  = (M)  ' e e (; 
z z z y 0 0 0 w 

where  the  prime  (’)  symbol  indicates  quantities  nondimensionalized  on  the 
basis  of  U,  b,  and  p/2,  and 

z',z',z'  * vertical  heave  acceleration,  velocity,  and  displacement, 

respectively,  positive  down 
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pitch  angular  acceleration,  velocity,  and  dia- 
placement,  respectively,  positive  bow  up 

mass  of  the  boat 

pitch  moment  of  inertia  of  the  boat  about  the  y-axis 
through  the  center  of  gravity 

wave  encounter  frequency  kb  ^1  + 


(Z)  ' , (M)  ' = complex  amplitudes  of  wave  excitation  force  and 

wave  excitation  moment,  respectively,  about  the 
center  of  gravity 

b = boat  beam 

k = wave  number,  2n/L 

w 

c = wave  celerity 

= wave  length 

U = steady  reference  boat  speed. 

The  coefficients  of  the  motion  variables  such  as  Z..',Z.',  etc.,  are  the 

Z 2 

stability  derivatives  of  the  boat.  There  are  12  of  these  in  the  previous 
equations  and  24  when  the  surge  degree  of  freedom  is  included.  The 
equations  for  these  coefficients  were  derived  in  Reference  4.  The 
equations  for  the  wave  excitation  force  and  moment  are  derived  in 
Appendix  A by  treating  the  craft  as  a slender  body  with  a correction  factor 
for  three-dimensional  effects.  These  equations  are  conveniently  written 
in  terms  of  five  components.  From  Equations  (120)  and  (122)  of 
Appendix  A,  they  are 


<ZV  e 


-iu)  't' 
e 


- - [<FlV  + <F2>v'  + <F5)w'  + <V»’  + (Fh>„’]  <=“  T »> 


-iui  't' 

«V  - e ■ («!)„'  * «2V  + («5V  * <Vv’  ♦ <Vv' 


The  definitions  of  the  individual  components  are  given  as  follows,  and  the 
appropriate  equations  for  regular  head  and  following  waves  are  found  in 
Appendix  A. 
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<FiV‘(MiV 

(F2V’<M2V 


(F,)  ',(M.)  ' 
5 w 5 w 


<Vw''<Vw' 


<rh>»'-<Mh>»' 


Integrated  effect  of  the  wave  slope  on  the  boat  to 
chine  lanersion;  see  Equations  (74)  and  (102) 

integrated  effect  over  the  length  of  the  boat  of 
the  wave  orbital  velocity;  see  Equations  (82) 
and  (103) 

integrated  effect  over  the  length  of  the  boat  of 
the  wave  orbital  acceleration;  see  Equations  (87) 
and  (104) 

integrated  effect  over  the  length  of  the  boat  of 
the  ambient  perturbation  pressure  due  to  the  waves; 
see  Equations  (111)  and  (115) 

effect  of  change  in  boat  wetted  area  due  to  wave 
elevation;  see  Equations  (118)  and  (119). 


All  of  the  previously  defined  quantities  are  derived  for  the  boat 
moving  unperturbed  at  constant  speed  U in  regular  head  or  following  waves. 
The  effect  of  other  wave  headings  is  readily  obtained  as  indicated  in 
Appendix  A.  The  hydrodynamic  and  inertial  forces  and  moments  resulting 
from  the  boat  motions  are  given  by  the  left-hand  side  of  Equations  (1) 
and  (2). 

NATURAL  FREQUENCIES  AND  DAMPING 

By  setting  the  right-hand  side  of  Equations  (1)  and  (2)  equal  to  zero 
and  solving  for  the  boat  response  to  an  initial  perturbation,  we  obtain  a 
solution  which  describes  the  dynamic  characteristics  of  the  boat  in  calm 
water. 


z'(t) 


0(t)  * 0^  e 


' o'  t' 

+ z2  e + 

o 't' 

+ ©2  e + . 


where  z^,  z2,....0^»  02»..«  are  constants  which  depend  on  the  initial  con- 
ditions. The  o'  terms  determine  the  character  of  the  time  history  response 
of  the  boat  to  any  small  disturbance.  Four  values  of  a'  are  obtained  from 
the  roots  of  the  resulting  characteristic  equation. 
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The  roots  of  these  equations  may  be  real  or  complex  conjugate  pairs. 

In  either  case  it  is  seen  from  Equations  (5)  and  (6)  that  if  any  root  has 
a positive  real  part,  the  transient  response  increases  without  limit,  and 
the  boat  is  considered  unstable  in  the  linear  sense. 

In  general  a complex  pair  of  roots  represents  an  oscillatory  mode; 
e.g.,  for  the  root  pair  o'  * o^'  + i Oj',  the  z'  response  is 

z'  * e (cx  cos  oI't'  + c2  sin  o 't')  (9) 

where  c^  and  c^  are  real  constants  which  are  determined  by  the  initial 
conditions.  The  magnitude  of  the  imaginary  part  of  the  root  cr^'  is  the 
nondimensional  natural  frequency  of  the  modal  motion.  In  dimensional  form, 
the  natural  frequency  and  period  are 

0 j = Oj*  g rad/sec  (10) 
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The  effect  of  the  real  part  of  the  root  O^'  may  be  illuatrated  by  computing 
the  time  for  a tranaient  disturbance  to  either  halve  or  double  Itself  in 
magnitude.  Thus,  if  oR'  is  negative,  the  envelope  of  the  disturbance  will 
be  halved  when 
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It  follows  that  the  time  for  the  disturbance  motion  of  each  mode  to  halve 
or  double  itself  is 
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0.69/o^  sec 


Another  useful  measure  of  damping  of  oscillatory  modes  is  the  damping  ratio 
which  is  directly  related  to  the  rate  of  decay  of  disturbance 
oscillations.  It  is  given  by 
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In  the  vicinity  of  the  resonant  encounter  frequency,  in  waves,  the  damping 
ratio  is  also  inversely  related  to  the  amplification  ratio  of  the  boat 
response.  Values  of  £ between  0.6  and  1.0  are  usually  considered  to  give 
well-damped  modes.  Values  less  than  about  0.4  are  generally  considered  to 
produce  underdamped  modes.  Although  the  foregoing  may  provide  a rough 
indication  of  the  vertical  plane  dynamic  characteristics  of  the  boat,  a 
dynamic  motions  analysis  is  required  for  any  detailed  study. 
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STEADY  STATE  MOTION  IN  WAVES 

The  steady  state  solution  of  Equations  (1)  and  (2)  is  a simple 
harmonic  motion  in  heave  and  pitch  and  has  the  following  form 
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where  z^'  and  0g  are  the  complex  amplitudes  of  the  heave  and  pitch  motion. 
Substituting  these  equations  into  Equations  (1)  and  (2)  leads  to  the 
following  equations  for  z^'  and  6q. 
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The  solution  to  these  equations  is 

<Z)w’  B2  ‘ (M)~'  A< 


'w  2 


A1  B2  " A2  B1 


- V + i V 


(19) 


(M)  ’ A.  - (Z)  1 B. 

ft  . W 1 V 1 _ n . . n 

60  Ax  B2  ~ A2  B1  • 0R  + 1 0! 


(20) 


Alternately,  in  terms  of  the  amplitude  and  phase  angle 
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The  denominator  of  Equations  (19)  and  (20)  is  the  characteristic 
quartic  (Equation  7)  where  iw  ' replaces  a1.  This  may  be  written  in  terms 
of  the  stability  roots.  When  this  is  done.  Equations  (19)  dnd  (20)  become 
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where  A is  given  by  Equation  (8)  and 


Q - (iu>e'  - a1’)  (iu>e'  - 02’)  (iu)e'  - o3')  (ioo^ * - a 4') 


For  the  planing  hulls  investigated  there  was  always  one  dominant  pair  of 
complex  conjugate  stability  roots.  If  the  damping  ratio  of  the  mode  is 
low,  one  may  expect  large  motions  when  the  wave  encounter  frequency  is 
near  the  modal  natural  frequency.  This  is  clear  from  the  following  form  of 
the  expression  for  modulus  of  Q in  Equations  (23)  and  (24). 
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where  u>m'2  - Oj 


oscillatory  mode,  and  t ^ 


is  the  square  of  the  undamped  natural  frequency  of  the 
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is  the  damping  ratio  of  the 


Similar  definitions  apply  to  the  o',  o'  modes. 


oscillatory  mode  (£,  < 1). 
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except  in  most  cases  ^ > 80  that  the  second  mode  is  not  oscillatory. 

It  is  seen  that  the  first  term  can  have  a sharp  minimum  at  the  resouant 
value  of  u)e' , depending  on  how  small  is.  If  the  remaining  terms  and 
the  wave  force  vary  slowly  with  frequency  in  the  vicinity  of  this  minimum, 
then  the  heave  and  pitch  amplitude  might  be  expected  to  have  a well- 
defined  resonance  peak  near  this  frequency.  The  encounter  frequency  at 
which  this  would  occur  is  given  by 
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where  u>DN'  = ^ |o^'  - o ^ | is  the  damped  natural  frequency  of  the  mode. 

It  is  seen  that  the  value  of  this  so-called  resonance  frequency  is  smaller 
than  both  the  damped  and  undamped  natural  frequencies.  If  Equation  (26) 
is  substituted  for  in  Equation  (25),  the  magnitude  of  f Q j is  approxi- 
mately proportional  to  for  small  and  the  response  z^'  and  9^  at 
resonance  is  roughly  inversely  proportional  to  Thus,  for  the 

previously  stated  conditions 
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COMPARISON  WITH  EXPERIMENT 

Since  Che  only  systematic  experimental  data  available  about  motions 

of  prismatic  planing  hulls  (Figure  1)  in  regular  waves  are  those  obtained 
2 

by  Fridsma,  the  theory  will  be  compared  with  these  data  only.  This  is 
not  particularly  restrictive  since  experiments  have  been  performed  with 
16  configurations  built  around  a rather  large  series  of  constant  deadrlse 
models  (Figure  1),  and,  as  stated  earlier,  measurements  were  obtained  for 
effects  of  deadrise,  trim,  loading,  length-to-beam  ratio  and  speed  on  the 
pitch- and- heave  amplitude  and  the  phase  angle  in  regular  head  waves.  To 
facilitate  comparisons  of  the  theory  with  these  measurements,  the  ampli- 
tudes of  the  heave  and  pitch  motion  are  nondimens ionalized  by  the 
amplitudes  of  the  wave  height  and  wave  slope,  respectively.  These  so- 
called  response-amplitude  operators  (RAO)  are  then 
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heave  RAO 


- pitch  RAO  (29) 

Since  the  sense  of  positive  boat  displacement  was  assumed  opposite  to  that 
for  wave  surface  elevation  in  the  equations  of  motion  (Appendix  A),  the 
phase  lag  e^,  when  both  displacements  are  taken  in  the  same  sense,  is 
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The  phase  lead  in  pitch  is  given  by 


ee  “ 2lT  " ae 


The  previously  described  quantities,  Equations  (28)  to  (31),  are  plotted 
against  a modified  version  of  the  nondimensional  wave  number 
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The  factor 


was  determined  empirically  so  as  to  minimize  the 


effect  of  the  boat  length- to-beam  ratio,  £/b  (or  rather  the  radius-of- 
gyration — to — beam  ratio)  and  the  load  factor  (or  the  mass)  on  the  value 
of  at  the  resonant  encounter  frequency. 

The  principal  geometric  characteristics  of  each  of  the  16  configu- 
rations investigated  are  given  in  Table  1.  The  values  of  the  non- 
dimensional  longitudinal  position  of  the  center  of  gravity  from  the 
transom  (LCG/b)^  shown  in  the  table  were  obtained  from  interpolations  of 

cross  plots  of  LCG  versus  trim  angle  curves  obtained  from  smooth  water 

2 

experiments.  These  were  estimated  by  Fridsma  to  give  the  running  trim 
angles  T_  shown  in  the  table. 

For  the  theoretical  predictions,  it  was  necessary  to  calculate  the 
running  trim.  This  was  determined  for  all  configurations  from  the 
previously  mentioned  values  of  (LCG/b)^.  Although  most  of  the  calculated 
trim  angles  T£  were  within  approximately  one-half  degree  of  the  desired 
t£  values,  a few  were  about  1 degree  lower.  Since  the  trim  angle  plays 
the  major  role  in  determining  the  dynamic  characteristics  of  the  boat,  it 
has  been  necessary,  in  these  cases,  to  make  calculations  for  a second  LCG 
position  selected  to  give  a calculated  trim  angle  closer  to  the  desired 
value.  The  resulting  determinations  are  shown  in  Table  IB  and  designated 
by  subscript  2. 

For  each  configuration  theoretical  calculations  were  made  of  the 
pitch  and  heave  RAO's  and  the  corresponding  phase  angles  as  defined  by 
Equations  (28)  to  (31)  for  wavelengths  ranging  between  1 and  60  boat 
lengths.  In  addition  the  stability  roots  for  each  configuration  were 
also  determined.  It  was  found  that  for  speed-to-length  ratios  of  4 and 
6 there  was  only  one  pair  of  complex  conjugate  roots  in  virtually  every 
case.  For  a speed-to-length  ratio  of  2,  however,  there  were  two  pairs 
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of  complex  conjugate  roots  in  every  case.  This  presumably  was  due  to  the 
importance  of  buoyancy  forces  at  these  lower  speeds.  Calculations  were 


also  made  of  the  value  of  C 
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This  is  the  value  of  C,  for  which  the 


encounter  frequency  is  equal  to  the  so-called  resonance  frequency  w ' as 

K 


determined  from  Equation  (26) 

Equation  (32)  is  given  by  the  following  equation 


This  occurs  for  head  waves  when  L /b  in 
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where  C = --- - is  the  speed  coefficient. 
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These,  together  with  the  estimated  damping  ratios,  are  shown  in  Table  1. 


It  is  seen  that  the  values  of  C^R  for  nearly  all  cases  lie  between  0.075 


and  0.10.  The  effect  of  increased  trim  angle,  decreased  load  coefficient, 
and  increased  speed  are  seen  to  decrease  the  damping  ratios.  The  effect 
of  length-to-beam  ratio  does  not  appear  to  be  significant.  The  effect  of 
deadrise  appears  to  be  small  at  a speed- to-length  ratio  of  4;  however,  at 
a speed- to- length  ratio  of  6,  the  damping  ratio  decreases  with  decreasing 
deadrise  angle  for  similar  running-trim  angles.  The  10-degree-deadrise 
configuration  J,  having  an  estimated  damping  ratio  of  0.06  at  a trim  angle 
just  under  4 degrees,  is  seen  to  have  the  lowest  damping  ratio  by  far.  It 
will  be  seen  that  this  condition  resulted  in  the  most  violent  motion  of  all. 

The  results  of  the  theoretically  determined  RAO's  and  phase  angles 
are  plotted  in  Figures  2 through  17.  The  figures  for  each  speed-to-length 
ratio  are  discussed  separately. 


SPEED-TO-LENGTH  RATIO  EQUALS  4 

Most  of  the  parametric  study  has  been  performed  at  the  speed-to-length 
ratio  of  4 because  this  speed  is  more  typical  in  planing  craft  operation 
than  speed-to-length  ratios  of  2 and  6.  Figures  2 to  9 show  comparisons 
of  the  measured  and  computed  values  of  the  pitch  and  heave  RAO  and  phase 
angle  as  a function  of  the  modified  nondimensional  wave  number  parameter  C^. 
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Also  sho%m  on  Che  abscissa  Is  a scale  of  wavelength- to-beam  ratio. 

The  solid  curves  represent  RAO's  and  phase  angles  as  calculated  on  the 


basis  of  the  model  (LCG/b)^  values  shown  in  Table  1A.  In  some  cases  broken 


curves  are  also  shown.  These  represent  the  effect  of  a small  change  in  the 
LCG  position  and  trim  angle  as  listed  in  Table  IB.  The  data  points  are 
indicated  by  the  symbols  shown.  The  experimental  wave  height- to-beam  ratio 
was  0.111  for  all  configurations,  except  the  one  shown  in  Figures  2a  and 
2b,  where  the  effect  of  wave  height  was  also  investigated.  It  may  be  seen 
from  this  latter  figure  that  the  effect  of  wave  height  does  not  appear  to 


be  large,  except  perhaps  at  C,  - 0.15.  However,  it  will  be  seen  later 


that  this  is  not  necessarily  true  for  all  configurations. 

Examinations  of  Figures  2 through  9 shows  that  the  motions  reach  their 
maximum  amplitudes  very  close  to  the  predicted  resonant  encounter  frequency, 
corresponding  to  the  value  of  CAR  given  in  Table  1A.  It  is  also  seen  that 
the  slope  of  the  phase  angle  curve  tends  to  be  a maximum  near  this  value 
of  C^,  as  one  might  expect.  The  overall  agreement  between  the  calculated 
and  measured  RAO's  and  phase  angles  is  seen  to  be  quite  good  even  for  the 
cases  in  which  the  conditions  of  the  experiment  were  rot  simulated  by  the 
theory  in  all  respects.  For  example,  the  running- trim  angle  for  Configu- 
ration E was  estimated  to  be  4.8  degrees,  while  the  model  value  was 
reported  to  be  6.0  degrees.  A recalculation  using  an  LCG  position  of  1.50 
beams  instead  of  the  original  1.73  beams  gave  a calculated  running- trim 
angle  of  5.8  degrees  and  better  agreement  with  the  data,  as  shown  by  the 
broken  curves  in  Figures  3a  and  3b.  Although  the  assumed  LCG  position  was 
further  aft  than  on  the  model,  this  discrepancy  had  much  less  influence  on 
the  boat  dynamics  than  the  running-trim  angle.  It  is  of  interest  to  note 
from  Table  IB  that  the  increased  trim  angle  resulted  in  a decrease  in  the 
damping  ratio  from  0.39  to  0.29,  thus  accounting  for  the  sharper  peak 
response.  It  is  also  worth  noting  that  the  phase  angles  in  Figures  3a  and 
3b  were  unaffected  by  the  trim  angle  change  to  a value  of  C^  of  0.15. 
However,  for  larger  values  there  is  a larger  phase  lag,  though  the  motion 


amplitudes  are  not  greatly  affected.  Fridsma  made  the  same  observations 
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on  the  basis  of  comparing  the  data  from  Configurations  A with  those  from  E; 
see  Figures  2 and  3.  These  configurations  have  estimated  damping  ratios 
of  0.47  and  0.29,  respectively.  However,  when  a boat  is  in  a more  stable 
configuration  due,  for  example,  to  increased  loading  or  smaller  running 
trim  angle,  the  damping  ratio  is  higher,  and  the  effect  of  a small  increase 
in  trim  angle  is  significantly  smaller.  This  is  shown  in  Figures  4,  7, 
and  9. 

In  a few  cases,  where  part  of  the  bow  was  estimated  to  be  in  water,  the 
agreement  with  theory  is  better  than  one  might  expect,  since  the  present 
theory  does  not  include  this  effect.  This  was  especially  true  for 
Configurations  K and  P;  see  Figures  6 and  9.  The  former  had  a 30-degree- 
deadrise  angle  and  the  latter  a length-to-beam  ratio  of  4. 

SPEED-TO-LENGTH  RATIO  EQUALS  6 

At  a speed-to-length  ratio  of  6,  the  planing  craft  becomes  less  stable 
and  the  motions  more  severe.  According  to  Table  IB,  the  damping  ratios  at 
this  speed  are  the  lowest.  It  is  seen  that  the  damping  ratios  of 
Configurations  B,  J,  and  M decrease  with  a decrease  in  deadrise  for 
running- trim  angles  slightly  under  4.0  degrees.  On  the  other  hand. 
Configuration  G,  because  of  its  higher  loading  has  the  highest  damping 
ratio  even  though  the  trim  angle  is  close  to  5.0  degrees. 

The  effect  of  wave  height  was  investigated  with  Configuration  B 
(Figure  10)  only.  It  is  seen  from  this  figure  that  the  nonlinear  effects 
are  more  pronounced  than  at  the  lower  speed.  This  is  especially  true 
near  the  resonant  encounter  frequency  and  may,  in  part  at  least,  account 
for  the  measured  peak  RAO's  being  lower  than  the  theoretical  ones. 

Reasonably  good  agreement  was  obtained  between  the  computed  and  measured 
RAO's  for  the  more  highly  loaded  20-degree-deadr ise  boat;  see  Figure  11. 

The  predictions  for  the  30-degree-deadrise  boat  (Figure  12)  were  not  as 
good,  particularly  for  the  pitch  RAO.  This  may,  in  part,  be  due  to  the 
fact  that  a portion  of  the  model  bow  was  in  the  water.  The  poorest 
agreement  with  the  predictions  of  RAO  were  for  the  10-degree-deadrise  boat; 
see  Figure  13.  The  possible  reasons  for  this  will  be  discussed  later. 
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For  every  data  point,  except  one,  by  using  the  theory  the  variation 
of  phase  angle  with  encounter  frequency  is  predicted  quite  well.  This  is 
especially  evident  from  the  solid  phase  angle  curves  shown  in  Figures  10 
through  13,  which  were  computed  from  the  original  LCG  position.  The 
reason  for  the  discrepancy  at  ■ 0.10  in  Figure  13  is  not  apparent. 
Calculations  for  the  more  aft  position  of  LCG  generally  show  an  increase  in 
the  phase  lag  at  C^  > 0.15  which  is  larger  than  that  at  a speed-to- length 
ratio  of  4.  This  indicates  that  the  boat  response  is  more  sensitive  at 
this  speed  than  at  the  lower  speed  to  the  shift  in  LCG  position  required 
to  correct  the  estimated  running  trim.  It  appears  that  an  Improved 
method  of  estimating  running-trim  angles,  within  the  context  of  the  present 
theory,  would  be  desirable  and  would  contribute  toward  removal  of  this 
type  of  discrepancy,  not  only  in  the  phase  angle  but  also  in  RAO 
predictions. 

As  noted  earlier,  poorest  agreement  between  theory  and  data  was  ob- 
tained for  RAO's  of  Configuration  J;  see  Figure  13.  According  to  Table  IB 
this  boat  had  an  estimated  damping  ratio  of  approximately  0.06  and  was 
therefore  not  far  from  a porpoising  condition.  This  is  indicated  also  by 
the  sharply  tuned  peak  in  measured  RAO's  and  the  steep  slope  in  the  phase 
angle  data  curve.  Although  the  theoretical  phase  angle  curve  is  shifted 


to  the  right,  it  is  apparent  that  it  exhibits  a variation  with  C,  in  the 


vicinity  of  resonance  similar  to  the  data.  This  suggests  that  the 
theoretical  prediction  of  the  damping  ratio  is  reasonable.  Since  the  peak 
RAO  is  roughly  inversely  proportional  to  the  damping  ratio  (Equation  (27)) 
a small  error  in  damping  ratio,  when  this  quantity  is  less  than  0.1,  can 
lead  to  large  discrepancies  in  RAO.  It  also  becomes  more  difficult  to 
find  the  resonant  peak  experimentally.  One  would  expect  to  find  it  close 
to  the  region  of  maximum  slope  in  the  phase  angle  curve.  Since  there  was 
no  recorded  data  point  in  the  interval  of  wavelengths  between  15  and 
20  beam  widths,  where  the  phase  angle  changed  by  more  than  100  degrees,  it 
is  possible  that  the  resonant  condition  was  not  measured.  It  is, 
therefore,  also  possible  that  peak  RAO's  are  actually  higher  than  those 
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shown  by  data  points  in  the  figures.  Even  if  this  were  so  it  would  be 
expected  that  nonlinear  effects  would  limit  the  response  to  something  less 
than  predicted. 

It  is  also  seen  in  Figure  13  that  a second  peak  in  RAO's  occurred  for 
a wavelength  of  approximately  7.5  boat  widths.  According  to  Reference  2 
the  model  was  reported  to  rebound  from  a wave  crest,  completely  "fly  over" 
a second  wave  crest,  and  land  on  the  third.  It  thus  had  an  actual 
excitation  frequency  of  one-half  that  normally  expected.  This  is  very 
close  to  the  resonant  frequency  of  the  craft,  which  presumably  occurs  at 
a wavelength  between  15  and  20  beam  widths.  Although  the  phenomenon  is 
undoubtedly  due  to  nonlinear  effects,  it  is  doubtful  that  it  would  have 
occurred  if  the  boat  were  not  in  a highly  tuned  condition  due  to  the  fact 
that  it  was  not  far  from  a porpoising  condition.  It  should  be  remarked, 
however,  that  it  is  somewhat  surprising  that  the  resonant  frequency  of  the 
boat  seems  to  have  been  altered  very  little,  although  the  change  in  the 
flow  condition  on  the  boat  is  apparently  very  large. 


SPEED-TO-LENGTH  RATIO  EQUALS  2 

The  speed-to-length  ratio  of  2 represents  the  prehump  condition  where 
buoyancy  forces  are  predominant,  and  significant  side  wetting  above  the 
chines  occurs.  Furthermore,  at  this  speed  the  hydrodynamic  coefficients 
may  be  Froude  number  and  frequency  dependent.  Since  the  theoretical  pre- 
dictions of  the  hydrodynamic  forces  assume  no  side  wetting  or  Froude 
number  dependence,  it  should  be  expected  that  the  motion  predictions  would 
suffer  somewhat  because  of  this.  Furthermore,  it  was  estimated  that  most 
of  the  bow  of  each  of  the  boat  models  tested  at  this  speed  was  in  the 
water.  As  noted  earlier,  the  present  theory  does  not  include  the  effect 
of  bow  immersion.  Despite  these  shortcomings,  the  theory  was  exercised, 
and  comparisons  were  made  with  the  experimental  data  for  all  of  the 
configurations  tested. 


As  mentioned  before,  all  of  these  configurations  had  two  pairs  of 
complex  conjugate  roots,  representing  two  natural  frequencies.  However, 
one  of  them  was  a low-frequency,  highly  damped  mode.  The  values  of 
and  € for  the  other  mode  are  shown  in  Table  1.  It  is  seen  that  all  of 
the  configurations  (C,  D,  H,  and  L)  are  predicted  to  have  a fair  amount  of 
damping.  Comparisons  of  the  theoretical  predictions  of  RAO  and  phase 
angle  with  the  experimental  data  for  the  10-  and  20-degree-deadrlse  boats; 
i.e..  Configurations  C,  D,  and  H,  shown  in  Figures  14  through  16,  are 
better  than  might  have  been  expected.  The  predictions  for  the  30-degree- 
deadrise  model  appear  to  be  worse  than  the  others.  This  is  probably  due 
to  the  fact  that  this  boat  ran  with  the  greatest  amount  of  bow  immersion. 


CONCLUSIONS 

The  theoretical  method  proposed  in  the  present  report  appears  to 
provide  reasonably  good  predictions  of  the  linear  response  characteristics 
of  high-speed  planing  boats  in  waves.  Accuracy  of  the  predictions  of 
RAO's  and  phase  angles  is  best  around  a speed-to-length  ratio  of  4.  This 
is  fortunate  since  this  is  the  most  typical  speed  for  planing  craft 
operation.  Although  the  predictions  of  the  phase  angles  at  a speed-to- 
length  ratio  of  6 are  good,  there  is  a tendency  to  overestimate  resonant 
RAO's  in  most  cases.  This  problem  appears  to  be  most  severe  when  the 
boat  is  operating  near  its  porpoising  condition.  Part  of  the  reason  for 
overestimation  may  be  due  to  the  fact  that  nonlinear  effects  are  not 
accounted  for  in  the  theory.  In  addition  it  is  felt  that  a more  precise 
method  for  estimating  the  running- trim  angle,  in  the  context  of  the 
present  theory,  would  lead  to  a further  improvement  in  these  predictions. 

Since  the  present  theory  was  derived  for  prismatic  planing  hulls,  it 
does  not  take  into  account  detailed  variations  in  hull  geometry. 
Furthermore,  it  is  a linear  theory,  and  a number  of  simplifications  and 
approximations  have  been  incorporated  in  it.  It  is  anticipated  that 
further  development  to  eliminate  these  shortcomings  will  take  place. 
Nevertheless  it  is  considered  to  be  a useful  design  tool  in  its  present 
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APPENDIX  A 

FORCE  AND  MOMENT  DUE  TO  UAVES  AND  THE  EQUATIONS  OF  MOTION 
FORCE  AND  MOMENT  DUE  TO  UAVES 

Ic  is  assumed  that  the  craft  is  proceeding  in  straight-line  motion 
at  constant  speed  U and  arbitrary  heading  angle  f?h  through  regular  waves . 
Perturbations  from  the  straight-line  motion  are  assumed  to  be  zero  and 
only  the  force  and  pitch  moment  excitation  in  the  vertical  plane,  due  to 
waves,  are  to  be  determined  in  this  section.  Since  the  theory  is  con- 
cerned mainly  with  the  high-speed,  low-aspect  ratio  condition,  it  is 
assumed  that  the  craft  may  be  treated  as  a slender  body  with  a three- 
dimensional  flow,  or  aspect-ratio  correction,  and  that  unsteady  effects 
are  small. 

The  total  hydrodynamic  force  on  the  craft  is  taken  as  the  sum  of 
the  time  rate  of  change  of  the  transverse  momentum  in  the  vertical  plane 
imparted  to  the  water  by  the  presence  of  the  moving  craft  and  the  ambient 
perturbation  pressure  due  to  the  waves.  Determination  of  the  first  of 
these,  the  dynamic  part  of  the  force,  initially  follows  the  lines  used  by 

A 

the  author  in  deriving  the  equations  for  porpoising  planing  craft. 

This  type  of  analysis  was  originally  suggested  in  1924  by  Munk^  and 
Jones6  in  connection  with  the  analysis  of  airships  and  slender  wings  and 
has  more  recently  been  generalized  by  Bryson7  for  slender  finned  missiles. 


3Munk,  M.M. , "The  Aerodynamic  Forces  on  Airship  Hulls,"  National 
Advisory  Committee  for  Aeronautics  Report  184  (1924) . 

6Jone8,  R.T.,  "Properties  of  Low-Aspect-Ratio  Wings  at  Speeds  Below  and 
Above  the  Speed  of  Sound,"  National  Advisory  Committee  for  Aeronautics 
Report  835  (1946). 

7Bryson,  A.E.,  Jr.,  "Stability  Derivatives  for  a Slender  Missile  with 
Application  to  Wing-Body-Vertical  Tail  Configuration,"  Journal  of 
Aeronautical  Sciences,  Vol.  20,  No.  5,  pp.  297-308  (1953). 


21 


It  has  been  applied  to  the  problem  of  the  translational  impact  of 


8 9 10 

seaplanes  on  water  by  Mayo  and  others.  * 


The  equations  for  the  wave  force  and  moment  are  described  with 
respect  to  orthogonal  axes  Ox,  Oy,  and  Oz  with  origin  fixed  at  the 
center  of  gravity  of  the  boat.  The  axis  Ox  is  in  the  fixed  direction  of 
the  forward  reference  speed  U of  the  boat.  The  axis  Oy  is  starboard, 
and  Oz  is  vertically  down;  see  Figures  18  and  19. 

The  flow  over  the  hull  is  assumed  to  occur  in  transverse  planes  which 
are  fixed  in  space  and  oriented  normal  to  the  keel.  The  momentum  of  each 
layer  of  water  transverse  to  the  keel  is  ds,  where  u is  the  two- 
dimensional  added  mass  of  the  section  of  the  hull  at  point  s,  interacting 
with  the  wave  motion  in  the  section  of  flow  plane  of  length  ds,  and  £ is 
the  component  of  the  orbital  velocity*  normal  to  the  keel  at  that  section. 
The  coordinate  s is  measured  aft  from  the  foremost  immersed  station  along 
the  keel.  The  dynamic  part  of  the  normal  force  on  the  section  ds  of  the 
hull  is  the  time  rate  of  change  of  the  momentum  imposed  on  the  layer  of 
water  ds  by  the  presence  of  the  hull 


dF  - ^ (uO  ds 


(34) 


Both  u and  £ will  in  general  be  functions  of  the  longitudinal  position 
coordinate  x of  the  hull  section  and  time  t.  Since  the  time  derivative 
is  in  the  fixed  transverse  plane,  it  must  reflect  the  changing  coordinate 
x of  the  transverse  plane  with  time.  Thus 


8. 


Mayo,  W.L.,  "Analysis  and  Modification  of  Theory  for  Impact  of 
Seaplanes  on  Water,"  National  Advisory  Committee  for  Aeronautics  Report  810 
(1945). 


Milwitzky,  B. , "A  Generalized  Theoretical  and  Experimental  Investigation 
of  the  Motions  and  Hydrodynamic  Loads  Experienced  by  V-Bottom  Seaplanes 
During  Step-Landing  Impacts,"  National  Advisory  Committee  for  Aeronautics 
TN  1516  (1948). 
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Schnitzer,  E.,  "Theory  and  Procedure  for  Determining  Loads  and 
Motions  in  Chine — Immersed  Hydrodynamic  Impacts  of  Prismatic  Bodies," 
National  Advisory  Committee  for  Aeronautics  Report  1152  (1953). 


* 

This  is  reasonably  constant  over  the  depth  of  planing  craft. 
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(35) 


« w fI  w 

It  " at  “ u a? 


The  normal  force  on  the  entire  hull  due  to  the  orbital  motion  la 
obtained  by  integrating  Equation  (34)  along  the  wetted  length  of  the  hull 


1^  and  multiplying  by  a correction  factor  $(X)  to  account  for  the  three- 


dimensionality  of  the  flow. 


FD  - <P(X) 


_ 2.  2.  _ 
[ J ' H ds  + | “ If  d*j 


(36) 


The  sectional  added  mass  at  any  transverse  section  depends  on  the  sectional 
geometry  of  the  boat  and  the  magnitude  of  According  to  Reference  4 
the  transverse  sectional  added  mass  distribution  for  a prismatic  hull  form 
may  be  estimated  from  the  following  equations  for  the  separate  contributions 


£1 


2 2 


f-  f (8)  C 


0 < s < s , 
- cl 


(37) 


P Tt  b 


(1  - sin  8)  s 


cl  < 8 < *k 


(38) 


2 Bb  (C  - Cc2) 


sc2  - 8 - \ 


(39) 


where  = wetted  length  along  keel;  see  Figure  19 


c2 


1.  - 1 


m 


mean  wetted  length  of  hull 


/T  tan  8 . 

8cl  " TtaiTf  b ln  calm  water 


(40) 
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--  • * 


(0.5  (0.57  + 0.0016)  (tan  0/ (2  tan  t)  - 0.0066)  - 0.03)b  (41) 


The  quantity  s^  is  the  value  of  s at  which  the  chine  becomes  "effectively" 
immersed.  The  term  B,  which  is  Bobyleff's  function,  is  a function  of 
deadrise  angle.  The  sectional  added  mass  p at  any  section  is  the  sum  of 
the  contributions  p at  that  section.  It  follows  from  Equations  (37) 
through  (39)  that  we  may  put 


dp  • 3p  • 3p 

dt  = C 5 3s  COt  T 


so  that  Equation  (36)  becomes 


Fd  * <fr(A) 


("I  12 1 ds  + / “ S ds 

L o o 


where 


P n C f (8)‘ 


0 < s < s , 
- cl 


§B» 


s - < s < l, 
c2  - - k 


£ 

& 


The  corresponding  moment  is  obtained  by  integrating  the  product  of  the 
stripwise  force  and  the  moment  arm  from  the  center  of  gravity  (a-s) . 


Mp  = <D(A) 


_ t 

1 

0 


(a-s)  i2  ds  + 


(a-s)  PC 


••  1 

C ds 


6 ■ deadrise  angle,  degrees 
4c2  - ; at  s - ,c2 

b » boat  beam 


where  a is  the  value  of  s at  the  transverse  flow  plane  through  the  center 
of  gravity;  i.e.,  a - 1^  - LCG,  and  LCG  is  ths  distance  ftr*®  transom  to 
center  of  gravity;  see  Figure  18. 

• • • 

Determination  of  £,  C and  £ 

For  the  condition  of  the  boat  moving  with  constant  speed  and  trim 
angle  T through  waves,  we  may  write 

C ■ s tan  T + £ (47) 


where  the  first  term  is  the  calm  water  value,  and  £ is  the  contribution 

w 


from  the  waves.  With  the  aid  of  Equation  (35)  we  have 


C * U sin  T + ( 


since,  according  to  Figure  19  we  see  that 


- cos  T 


Also  since  U sin  T is  a constant 


C - C. 


The  values  of  5 , i,  and  C are  next  determined  from  the  velocity 
potential  <}>w  for  the  wave  motion.  This  may  be  written  with  respect  to 
the  moving  boat  axes  as  follows,  where  x,  y,  and  z are  measured  from  a 
point  directly  above  the  center  of  gravity,  and  the  real  part  is  to  be 
taken. 


-kz  ik  [x  cos  0 + y sin  0.  + (U  cos  0,_  - c)  tl 
<j>  - ihc  e e h h h 

w 


where  * angle  of  wave  crests  from  boat  x-axls;  8h  * 0 and  ir  for 

following  and  head  seas  respectively;  see  Figure  18 

h ■ wave  amplitude,  1/2  wave  height 

c * wave  celerity 

k ■ 2u/L  ■ wave  number 
w 

L - wavelength 
w 

t * time 


In  the  following  we  will  consider  only  the  head  and  following  sea  cases. 


We  then  have 


4>  - ihc  e“kz  e+ik[x  + (U  t c>  fcl 

w 


In  Equation  (51)  and  the  following,  the  upper  sign  refers  to  head  seas; 
the  lower  sign,  to  following  seas.  The  encounter  frequency  co  is 
given  by 

co  = k (U  + c)  0 

e - 


The  wave  elevation  is 


n(x,t) 


(^) 


+i(kx  + oj  t) 
e 


since  c = g/k.  The  horizontal  and  vertical  components  of  the  orbital 
velocity  are  given,  respectively,  by 


8 <p  -kz  +i(kx  + to  t) 

W — 6 

v = - — r — = + (oh  e e 
w dx 


9 d>  -kz  +i  (kx  + <o  t) 

w - - -r-*  - iuih  e e e 

w dz 


where  the  circular  frequency  of  the  orbital  motion  is  given  by 


to  * kc 
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The  expression  for  2^  »*y  be  obtelned  as  the  sum  of  the  components  of  the 
orbital  velocity  normal  to  the  keel. 


C ■ - v sin  T - v cos  T 
w w w 


Inserting  Equations  (54)  and  (55)  leads  to 


. -kz  +i(kx  - T + u t) 

C ■ - luih  e e 8 

w 


The  expression  for  5W  is  readily  found  from  this  equation  and  Equation  (53) 


C ” n e 
w 


By  operating  on  Equation  (57)  with  Equation  (35),  we  obtain  the 
expression  for  the  acceleration  term.  Equation  (49) . 


» -kz  +i(kx  - T + to  t) 

V V 4,  0 

C “ C * - w h e e 
w 


Evaluation  of 

If  we  substitute  Equations  (45)  and  (47)  through  (49)  into 
Equation  (44)  we  obtain,  on  discarding  second  order  terms  in  the  wave 
disturbance. 


F « F + F + F + F + F 
D 1 r2  3 4 f 5 


where 


Fx  - 4>(X)  p tt  f(B)2  U2  sin2 


* 2 4>(A)  p 7T  f(6)  U sin 


s 

* ! 


(s  tan  T + C^)  ds  (61) 


s 

T tan  T J 


s C ds 
w 


u 

r,'~t 


l 


The  first  term  represents  the  effect  of  the  variation  in  added  mass 
to  the  point  of  effective  chine  immersion  for  constant  normal  speed 
component  U sin  T.  The  next  two  terms  give  the  effect  of  the  orbital 
velocity.  The  fourth  term  is  the  calm  water  normal  force  contribution  over 
that  portion  of  the  hull  from  s^  to  the  transom.  The  last  term  gives  the 
force  due  to  wave  orbital  acceleration. 

In  Reference  A the  quantity  s^,  which  is  the  value  of  s at  which 

the  chine  just  becomes  effectively  immersed,  is  evaluated  on  the  basis 

11  12 
of  Wagner's  wave  rise  theory  and  an  empirical  correction  by  Shuford  for 

the  calm  water  case.  In  that  derivation,  the  value  of  £ is  s tan  T and  sci 

is  defined  by  the  following  equation 


H 

;1 


■Cl 


(f(8)  scl  tan  t)2  = y-  (1  - sin  8) 


(66) 


To  make  use  of  this  result,  the  expression  for  £w  in  Equations  (58)  and  (61) 
is  approximated  by  the  following  equation  in  the  range  s = 0 to  s *=  sc^* 


^Wagner,  H.,  "The  Phenomena  of  Impact  and  Planing  on  Water,"  National 
Advisory  Committee  for  Aeronautics  Translation  1366,  ZAMM  Bd  12,  Heft  A, 
pp.  193-215  (Aug  1932). 

12 

Shuford,  C.L.,  Jr.,  "A  Theoretical  and  Experimental  Study  of  Planing 
; rfaces  Including  Effects  of  Cross  Section  and  Plan  Form,"  National  Advisory 
Committee  for  Aeronautics  Report  1355  (1957). 
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4 


* 


5 - 8 

w 


Making  use  of  Equation  (58)  in  Equation  (67)  and  substituting  this  for  £ in 


Equation  (61) , we  obtain 


where 


s tan  T + Cw  - s tan  (T  + t ) + ^(0) 


t . > /1(kVT>  r . ,l  * 

W X , /cos  T (_e  J e 


* horizontal  distance  from  center  of  gravity  to  foremost 
Immersed  station  on  keel 


x , = s , /cos  T 
cl  cl 


The  first  term  on  the  right-hand  side  of  Equation  (68)  represents  the 
combined  effect  of  trim  angle  and  wave  slope  on  the  force  in  the  interval 


S “ 0 to  8 


The  second  term  is  the  effect  of  wave  elevation  on 


the  force.  This  effect  will  be  evaluated  later.  On  substituting  the  first 
term  on  the  right  side  of  Equation  (68)  into  Equation  (61)  and  making  use 
of  Equation  (66)  with  T replaced  by  T + T , we  obtain  to  the  first  order  in 
the  wave  term 


b 2 

4>(A)  p 7t  — (1  - sin  8)  U (sin  T cos  T - x ) 
4 w 


The  first  term  in  this  equation  is  the  calm  water  normal  force  contribution. 

When  this  is  added  to  we  obtain  the  total  dynamic  part  of  the  steady  state 

normal  force  F on  the  hull.  We  nondimensionalize  by  dividing  through  by 
2 2 

1/2  p U b and  obtain 
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FDS  * 1 +"x  T “ 8in  sin  T cos  T + CD  C sin  T COS  T * COS  ® (71) 

where  we  have  made  the  following  substitutions  from  Reference  4. 

<72> 

( A ) B = C_  cos^  T cos  6 (73) 

u > c 


CL  = cross  flow  drag  coefficient;  see  Table  2 
D,c 

X = mean  wetted  length-to-beam  ratio 


TABLE  2 - CROSS  FLOW  DRAG  COEFFICIENT* 


Section  Shape 

•CD,c 

V-Bottom,  constant  deadrise 

.1.33 

V-Bottom,  horizontal  chine  flare.. 

. .1.33 

+ 0.0147  e° 

V-Bottom,  vertical  chine  strips... 

. .1.60 

+ 0.0147  6° 

* 

Reference  12. 

If  we  nondimensionalize  the  wave  term  in  Equation 

(70)  we  obtain 

(F. ) ’ - - G T 

1 w 1 + A w 

(7 

where 

G = | (1  - sin  6) 

(7 

Equations  (62),  (63),  and  (65)  are  expressed  in  terms  of  the  orbital 
motion  at  some  depth  z . For  planing  boats  there  is  generally  very  little 

cl 

difference  whether  we  take  z at  the  keel  or  at  the  free  surface. 

it-  ■ww’ier*-  -^hk 


m:  :■> 


Nevertheless  we  approximate  It  by  the  equation 


r 


' r 


a 


3 


*a  ■ Y s sin  T 


(76) 


where  y is  taken  as  the  fraction  of  the  depth  to  the  keel  at  which  the 
orbital  motion  is  considered  to  act.  Since  integrations  are  to  be  carried 
out  from  the  foremost  wetted  portion  of  the  keel  aft,  we  substitute  for 
x in  Equations  (57)  and  (59)  the  following 


x ■ 


COS  T 


(77) 


Substituting  Equations  (76)  and  (77)  ir.to  Equations  (57)  and  (59)  gives 

(78) 


£ » - iu)h  e 

w 


+!(<(/  ~ kQ  s) 


, 2.  +1(*  ' k0  s) 
t,  * - 0)  h e 
w 


(79) 


where 


i|/»kx^-T+u)et 


kQ  - (k/cos  T)  (1  + i y sin  T cos  t) 


(80) 

(81) 


On  substituting  Equation  (78)  into  Equations  (62)  and  (63),  carrying  out 
the  integration,  and  nondimensionalizing,  we  obtain  with  the  aid  of 
Equations  (66),  (72),  (73),  and  (75)  the  following  expression  for  the 
normal  force  due  to  orbital  velocity. 


+ii|> 


<F2V  " F2*  + F3*  “ " 1 4 (v)  * {H1  cos  T + H2  8ln  T>  (82) 


where 


-rx  °8  2 [*"  0 cl  u ? i Vci>  - l] 

^Ko  scr 


for  8 > 2°  (83) 


X 7T 

1 + X 4 


for  8 < 2 


Cn  o -5  / +i  kn  s , +ik„  l' \ 

»2  = + 1 Tk~b  C0S  T C0S  S V ' e~  V 


The  selection  of  a 2-degree-deadrise  angle  in  Equations  (83)  through  (84) 

was  to  avoid  the  computational  problem  of  a near  zero  in  the  term 
2 

(kQ  s ) at  a deadrise  angle  close  to  zero.  It  is  assumed  that 
Equation  (84),  which  is  strictly  for  zero  deadrise  angle,  is  suitable  for 
angles  to  at  least  2 degrees. 

The  force  due  to  orbital  acceleration  is  obtained  by  substituting 
Equations  (37),  (38),  (39),  and  (79)  into  Equation  (65).  Thus 


9 +i^  _ 9 9 

F5  — <f)(X)  w h e f(B)  tanZ 


g 


2 k0S  . 
s e ds 


+ P^b-  (1-sin  8) 


+i  k_s 

0 , . p Bb  tan 

e ds  + * ^ 


A 

1 J 


(s-sc2)  e 


If  we  carry  out  the  integrations,  make  the  substitutions  indicated  by 
Equations  (66),  (72),  and  (73)  and  nondimensionalize,  we  obtain  finally 


-*  (?)2  ($)  * [t  - 


Y2-  + h4  + tan  x (H5  - H2  Xc2)  (87) 


1 


where 


X ° 2cl 
1 + X .3 

(k0  ’cl' 


.-1  ko*cl 


[2-<k0  ‘cl)2  ?12<k0  *cl>]-2  }'“»>■  2° 


for  8 < 2 


H4  * - 1 1 + A 2 k„b  ( 


-1  k0  8cl  i1  k0 
e - e 


0 kj 
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In  suiranary  Che  dynamic  pare  of  Che  nondimens Iona 1 force  on  Che  boae 
moving  wich  consCanc  speed  and  Crim  Chrough  regular  head  or  following  waves 
is  given  by 

<V  ■ V + + <F2>«'  + (F5>w'  (51) 


The  effecC  of  Che  ambienc  pressure  and  increased  weCCed  area  due  co  Che 
waves  will  be  decermined  laCer. 

Eva lua Cion  of 

The  dynamic  pare  of  Che  momenc  abouc  Che  boaC  cenCer  of  gravicy  is 
obcained,  according  Co  EquaCion  (46),  by  mulciplying  each  inCegrand  in 
EquaCions  (61)  Chrough  (65)  by  Che  momenc  arm  (a  - s) . The  corresponding 
momenC  equaCion  becomes 


% - Mx  + M2  + M3  + M4  + M5 


where 


* 4>(X)  pn  f(8)2  U2  sin2  T tan  (T  + 
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M4  - 0(A)  B | U2  sin2 


Tb  J 


(a-s)  ds 


M5  ■ 4> (A)  I (a-s)  u ds 


On  carrying  out  the  integration  in  Equation  (93)  and  using  the  result 
of  Equation  (70)  one  readily  obtains 


M.  = 4>(A)  (1-sin  8)  U2  (sin  t cos  T - T ) ( 

w \ 


L - LCG  - -r  S 


3 cl 


) (98) 


The  last  term  in  parentheses  is  the  moment  arm  of  from  the  center  of 
gravity.  It  is  readily  evaluated  with  the  aid  of  Equations  (40)  and  (41) 


where  T is  replaced  by  T + T . Thus 


Z - LCG  - 4 s , 
k 3 cl 


tm  - LCG  - (o.  157  - Hirffrr ) b + b R(6)  (99) 


where  R(8)  = -0.0038  (0.57  + 0.0018)  - 0.0 3 
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Combining  this  with  Equation  (98)  and  discarding  terms  of  higher  order 
than  one  in  the  wave  components  yields 

Mx  ■ ♦ (*)  SL^~  <l-sin  8)  U2b  ^ain  T cos  T (Xfc  - Xg  - j Xcl ) 

( 

-(x-2(o.!57  - ^16)  ^1-Xg-R(6))  xj 


(100) 


where  X . is  s ,/b  for  calm  water.  The  first  term  in  the  brackets  is  the 
cl  cl 

calm  water  moment  contribution.  When  this  is  added  to  we  obtain  the 
total  dynamic  part  of  the  steady  state  moment.  We  nondlmenslonallze  this 
by  dividing  by  1/2  p U2b^  and  obtain 


V ■ TTT  0 sl"  T “•  T (*k  - xg  - ! Aci) 

+ C_  sin2  T cos2  T cos8  X ( 4 - X ) 
D,c  '2  g/ 


(101) 


where  Equations  (72),  (73),  and  (75)  have  been  made  use  of.  If  we  non- 
dlmenslonallze the  term  in  Equation  (100)  due  to  the  wave  slope  we  obtain 


(M.)  ' 
1 w 


1 + X 


2 (0.157) 


tan  8 


] T» 


(102) 


where  the  terms  in  8 were  dropped  since  they  were  found  to  have  very  little 
effect  on  the  accuracy  of  the  calculations. 

The  part  of  the  moment  (M2  + M^)  due  to  the  wave  orbital  velocity  is 
seen  to  be  expressible  as  the  sum  of  the  moment  of  the  force  (F.)  acting 

L W 

at  s * 0 and  the  contribution  obtained  from  the  integration  of  the  second 
term  in  the  integrands  of  Equations  (94)  and  (95) . On  carrying  out  this 
operation,  making  the  appropriate  substitutions  from  Equations  (66),  (72), 
(73),  and  (75)  and  nondimens ionalizlng,  one  obtains 


<m2)w’  - M2*  ♦ M3* 


‘f2V  <W  4 4 (f)  (t)  41"  3 '***  (1h5  - Tzn) 


That  part  of  Che  moment  My  due  to  Che  acceleration,  la  expressible  as 
Che  sum  of  three  Integrals  analogous  to  that  given  In  Equation  (86)  for 
Fy  On  carrying  out  the  integrations  Indicated  in  Equation  (97),  making 
the  appropriate  substitutions  as  before,  and  nondimensionalizing  one 


obtains 
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In  summary  Che  dynamic  part  of  Che  nondimenaional  moment  on  Che  boat 
moving  with  conatant  speed  and  trim  through  regular  head  or  following  waves 
is  given  by 

<V'  ■ "os'  + (MlV  + <M2>.’  * <MsV  (105> 


Normal  Force  and  Moment  Due  to 
Ambient  Pressure 


The  force  due  to  local  pressure  effects  Fp  is  obtained  by  integrating 


the  ambient  pressure  over  the  bottom  of  the  hull.  Thus 


dS 


(106) 


where  the  second  term  is  the  effect  of  the  wave  perturbation 


d <|> 


dt 


-kz  +i(kx  + a)  t) 
ghee 


(107) 


and  the  first  term  represents  the  static  buoyancy  force.  We  will  consider 


only  cases  for  speed  coefficient  Cy  at  more  than  0.5,  where  the  water 


breaks  clear  of  the  transom,  thus  ventilating  the  backside  of  the  boat  to 
the  atmosphere  and  removing  the  component  of  force  parallel  to  the  keel.  For 
this  case  the  force  due  to  ambient  pressure  may  be  assumed  to  act  normal  to 
the  keel. 

found  to  fit  existing  data  reasonably  well." 


The  following  expression  for  the  static  buoyancy  force  was 

13,14 


FBS  = K \ P 8 b3  A2  sin  T 


(108a) 


13. 


Hsu,  C.C.,  "On  the  Motions  of  High  Speed  Planing  Craft,"  Hydronautics 
Report  603-1  (May  1967). 


14, 


Brown,  P.W.,  "An  Experimental  and  Theoretical  Study  of  Planing  Surfaces 
with  Trim  Flaps,"  Davidson  Laboratory,  Stevens  Institute  of  Technology 
Report  1463  (Apr  1971). 
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This  is  equivalent  to  < p g b 


where  < is  an  empirical  correction  factor,  which  accounts  for  the  ventilation 

14 


effects.  A value  of  0.624,  obtained  from  data  on  prismatic  hulls,  has 
been  used  in  the  present  calculations  and  in  Reference  4.  Thus 
Equation  (108a)  becomes,  in  nondimens ional  form 


„ . 0.624  ,2  . 

fbs  " — T A sin  T 

CV 


(108b) 


For  the  sake  of  consistency  with  the  above  result,  the  same  value  for  < 
and  the  same  integration  limits  were  used  also  in  evaluating  the  force  due 


to  the  wave  perturbation  pressure  (F  ) . Thus  the  second  term  in 

r W 


Equation  (106)  becomes 


+1  Ip. 
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y . -(sin  T +i  cos  T)  k(s  - s _) 

(F)  - 0.624  pgbh  e d | e ds  (109) 

p W 


c2 


where  the  following  substitutions  have  been  made  in  Equations  (106) 
and  (107) 


dS  * b ds 


z = (s  - s^)  sin  T 


ip  . - kx , + u>  t 
d d e 


(110) 


*d  " *b  ' Sc2 


- s . COS  T 


X = - (s  - Sc2>  COS  T 


On  carrying  out  the  integration,  we  obtain  in  nondimensional  form 
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where  Cy  ■ U/4b 


Thus  the  total  normal  force  due  to  the  ambient  pressure  is 


V - fbs*  * <fpV 


(112) 


According  to  Reference  14  the  assumption  that  the  buoyancy  force  acts 
one-third  of  the  mean  wetted  length  forward  of  the  stern  gives  good 
agreement  with  data.  Accordingly,  the  nondimens ional  moment  about  the 
center  of  gravity,  due  to  the  buoyancy  force,  is 


V ’ TT  <x2  T>  (3  - Ag) 


(113) 


The  moment  due  to  the  wave  perturbation  pressure  is,  with  the  aid  of 
Equation  (109) 

i.. 


(Mp)w  * 0.624  p gbh  e 


+i 
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-(sin  T + i cos  T)  k(s-s  „) 
(s^-s)  e c ds 


c2 


(114) 


where  s.  ■ (A  - A )b 

d g 


Carrying  out  the  integration,  we  obtain  in  nondimensional  form 

<V„'  ■ <V.'  <x  - V + 5 1 


+ i(l|»d  + 2 T) 


k b 


(1  - H?)  (115) 


39 


a a min 


■p  W : 


where 


-k  X b(sin  T + i cos  t) 

* [1  + k X b(sin  T + i cos  T) } e 


Thus  Che  total  moment  about  the  center  of  gravity  due  to  the  ambient 
pressure  is 


<V  ‘ V * <Vw' 


(116) 


Normal  Force  and  Moment  Due  to  Change 
in  Wetted  Length 

It  is  recalled  that  only  part  of  the  force  and  moment  contributions 

arising  from  the  effect  of  wave  elevation  has  so  far  been  accounted  for. 

These  are  given  by  Equations  (74)  and  (102)  respectively,  and  represent  the 

effect  of  the  mean  wave  slope  T in  Equations  (67  through  69).  The  wave 

w 

elevation  term  £ (0)  of  Equation  (67),  which  must  also  be  included,  has  the 
w 

effect  of  increasing  the  wetted  length  of  the  boat.  The  effect  of  this  on 
the  force  may  be  determined  with  the  aid  of  Equations  (71)  and  (101),  which 
express  the  effect  of  the  mean  wetted  length  X and  trim  angle  T on  the 
steady  normal  force  and  moment  due  to  forward  speed  alone.  The  change  in 
mean  wetted  length  due  to  wave  elevation  will  be  approximated  from  the  free 
surface  elevation  at  the  point  x = x,. 
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A X 


+ i(k  x , + u)  t) 
, d e 

h e 

b sin  T 


(117) 


The  position  x^,  rather  than  ^x^  - , has  been  used  since  Xb  is  measured 

from  the  transom  to  the  point  s = s which,  we  have  seen  from 

c L 

Equation  (110),  is  directly  below  the  point  x = x^  on  the  free  surface. 

With  the  aid  of  Equations  (71),  (75),  and  (117)  we  may  estimate  the  non- 

dimensional  normal  force  (F,)  due  to  wave  elevation,  from  the 

n w 


► 


\ 
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following  expression 


x , 9 FDS*  . , 

Fh^w  “ 3 X A X 


which,  on  making  Che  appropriate  substitutions,  yields 


<FhV 
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(1  ♦ X)2  D’< 
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(118) 


The  nondimens ional  moment  about  the  center  of  gravity  due  to  wave 
elevation  is  obtained  in  an  analogous  manner  with  the  aid  of  Equation  (101) , 


9 “os' 

<Vw'  - TT-  A x 


(M,  ) ' = G sin  T cos  t =•  ( 

* W L (1  + A)2  \ 


X.  - X - ± X . 
k g 3 cl 


) +1M 


}+  i(k  + u>e  t) 

— 

b sin  T 


(119) 


The  next  step  is  to  incorporate  all  of  the  foregoing  results  for  wave 
excitation  into  equations  of  motion  of  the  boat. 


EQUATIONS  OF  MOTION 

The  equations  of  motion  of  planing  craft  in  the  longitudinal  degrees 
of  freedom  are  most  conveniently  expressed  in  terms  of  heave,  surge,  and 
pitch  motions.  Accordingly  the  normal  force  due  to  the  waves  are  resolved 
into  the  vertical  and  horizontal  components.  We  employ  the  convention  that 
the  nondimensional  surge  force  and  displacement  are  denoted  by  X'  and  x', 
respectively,  and  are  positive  in  the  direction  of  the  calm  water  speed; 
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Che  nondimens Iona 1 vertical  force  and  displacement  are  denoted  by  Z'  and 
z',  respectively,  and  are  positive  in  the  downward  direction  normal  to  the 
calm  water  free  surface.  The  nondimensional  pitching  moment  and  dis- 
placement angle  are  M*  and  0 and  are  positive  for  the  bow-up  direction  of 
rotation;  see  Figure  18. 

Since  the  motions  due  to  waves  will  be  determined  relative  to  the 

steady  state  attitude  and  draft  of  the  boat,  the  steady  state  hydrodynamic 

force  and  moment  make  no  contribution,  in  the  linear  sense,  to  the  boat 

motions.  This  is  discussed  further  in  Reference  4,  where  the  complete 

steady  state  equilibrium  equations  are  derived.  The  vertical  component 

of  the  linearized  wave  force  is  therefore  obtained  by  adding  the 

contributions  from  Equations  (91),  (112),  and  (118)  and  by  omitting  F ' 

D S 

and  fbs’- 

+i  uj  ' t ' 

(Z)w'  e 6 = " [<FlV  + (F2)w’  + (F5)w*  + (Vw’  + ^hV1  COS  1 (120) 

The  horizontal  component  of  the  wave  force  is 

+i  oi  ' t'  +i  a)  ' t’ 

(X)w'  e = (Z)w’  tan  T e e (121) 

The  moment  due  to  the  waves  is  obtained  by  summing  Equations  (105),  (116), 
and  (119)  and  by  omitting  M^'  and  M^'.  Thus 

+i  w ' t ' 

<M)„'  e £ " <Ml>w’  + <m2>w’  + (M5V  + (MP>w'  + (Vw’  (122) 

In  the  previous  equations  (X)  ' , (Z)  ',  and  (M)  ’ are  of  course  complex 

w w w 

constants. 

Equations  for  estimating  the  linearized  hydrodynamic  forces  on  the  boat 
due  to  perturbations  from  steady  state  motion  in  surge,  heave,  and  pitch 
have  been  derived  in  Reference  4 for  use  in  the  theoretical  prediction  of 
calm  water  porpoising.  The  same  equations  are  of  course  applicable  for 
calculating  the  linearized  hydrodynamic  force  and  moment  on  the  boat 
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arising  from  its  motion  response  to  wsves.  If  we  express  the  two  non- 
dimensional  force  components  and  the  nondlmensional  moment  due  to  the  boat 

motions,  as  linear  functions  of  the  nondimensional  perturbations  in  surge, 

4 

heave  and  pitch,  we  have 


(X')_.,  -X.  u + X u + X..  z + X.  £ + X 2 + Xx  0 + X;  0 + Xfl  9 (123) 

DM  UUZZZuuQ 

(Z')n  » Z.  u + Z u + Z..  z + Z;  z + Z z + ZK  0 + Z:  0 + Z„  0 (124) 

Bn  uuzzzono 

(M')nu  “ M.  u + M u + M»  z + M.  z + M z + Ms  0 + Ma  0 + Ma  6 (125) 

BM  UUZZZO0D 

where  u = x',  and  u * x'.  The  primes  used  to  denote  nondimensional  values 
have  been  omitted  on  the  right-hand  side  of  the  equations  for  convenience. 

The  coefficients  X^,  Xy,  etc.,  in  the  previous  equations,  are  the  (stability) 
derivatives  of  the  force  and  moment  with  respect  to  the  motion  perturbations 
u,  u,  etc.  The  dot  above  the  symbols  represent  derivative  with  respect  to 
time.  Theoretical  estimates  for  all  of  the  stability  derivatives  in  the 
previous  equations  are  derived  in  Reference  4. 

The  equations  of  motion  in  waves  may  now  be  written  in  terms  of 
Equations  (120  through  125)  as 


+1  w ' t ' 


»'*’  - <X'>bm  = <X>w’  e 6 

(surge  equation) 

(126) 

+ii)  ’t’ 

m*z'  - (Z’)BM=  (Z)w*  e e 

(heave  equation) 

(127) 

+i  w ’t' 

■y’8'  “ (M,)BM5  6 ^ 

(pitch  equation) 

(128) 

where  m'  = r is  the  nondimensional  boat  mass,  and  I ' = ~r  *■  r is  the 

± Pb3  * i Pb5 

nondimensional  pitch  moment  of  inertia  of  the  boat  about  the  y-axis  through 
the  center  of  gravity. 
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Equations  (126  through  128),  with  the  wave  force  and  moment  terms  set  equal 
to  zero,  are  the  equations  for  porpoising  stability.4  The  steady  state 
solution  of  Equations  (126  through  128)  is  a simple  harmonic  motion  in  the 
perturbations.  Thus 


+i  id  ' t ’ 

X'  = V e 6 


(129) 


+i  (d  *t* 
z’  = z0-  e e 


(130) 


9 = 0O  e 


+i  id  't' 
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(131) 


where  Xq',  zq' , 0Q  are  the  complex  amplitudes  of  the  displacement  and  pitch 
angle  perturbations.  If  we  substitute  Equations  (129  through  131)  into 
Equations  (126  through  128)  we  obtain  the  equations  for  the  amplitudes  and 
phase  angles  of  the  motions. 
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(133) 


(134) 


The  prime  on  the  brackets  indicates  that  each  term  is  primed.  The  solution 
to  these  equations  give  the  complex  amplitude  of  the  boat  motions 


v ■ v + 1 y • iy  i 


(135) 
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(137) 


where  the  subscripts  R and  I refer  to  the  real  and  imaginary  component. 
The  displacement  amplitudes  and  phase  angles  are  given  by  the  following 
equations 


xo’l  - 1 

/xr'2  + V2 

a * 

X 

tan  ^ 

(Xi’/V> 

zo’l  “ 1 

/Zr'2  + V2 

a ■ 
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(2IVV> 
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/ 6R2  + 6I2 

a0  “ 
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<vv 

Since  the  convention  for  heave  displacement  is  positive  down  while  the 
convention  used  for  the  change  in  free-surface  elevation  due  to  waves  is 
positive  up,  it  is  necessary  to  shift  the  phase  angle  for  the  heave  motion 
by  180  degrees  if  we  wish  to  determine  the  phase  angle  of  the  heave  motion 
with  respect  to  the  free-surface  displacement,  positive  up  in  both  cases. 

On  this  basis,  we  readily  find  that  the  phase  lag  of  the  heave  displacement 
with  respect  to  the  free-surface  elevation  n at  the  boat  center  of 
gravity  is 


e * a - n head  waves 
z z 


(138) 


e * it  - a 
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following  waves 


The  phase  leads  of  the  surge  amplitude  and  pitch  angle  with  respect  to  n 
at  the  same  point  are,  respectively 
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Figure  3 - Heave  and  Pitch  Responses  (V//f  ■ 4)  Configuration  E 
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Figure  4b  - Pitch  Response,  Speed-to-Length  Ratio  of  4 
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Figure  15  - Heave  and  Pitch  Responses  (V//£  - 2)  Configuration  D 


ENLARGED  DETAIL  OF  CHINE  IMMERSION 
(EFFECT  OF  WAVE  RISE  NOT  SHOWN) 


Figure  18  - Coordinate  System  and  Symbol  Illustrations 
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